GRAVITATIONAL WAVES AND DYNAMICS 
OF COMPACT BINARY SYSTEMS 



LUC BLANCHET 
Institut d'Astrophysique de Paris (C.N.R.S.), 
9^'" boulevard Arago, 75014 Paris, France 
E-mail: blanchet@iap.fr 

Part A of this article is devoted to the general investigation of the gravitational- 
wave emission by post-Newtonian sources. We show how the radiation field far 
from the source, as well as its near-zone inner gravitational field, can (in principle) 
be calculated in terms of the matter stress-energy tensor up to any order in the 
post-Newtonian expansion. Part B presents some recent applications to the prob- 
lems of the dynamics and gravitational- wave flux of compact binary systems. The 
precision reached in these developments corresponds to the third post-Newtonian 
approximation. 

1 Introduction 

The theory of gravitational radiation from isolated sources, in the context 
of general relativity, takes on importance because of the possibility of com- 
paring the theory with contemporary astrophysical observations. We have in 
mind the large-scale optical interferometers LIGO, VIRGO, GEO and TAMA, 
which should routinely observe the gravitational waves produced by massive 
and rapidly evolving systems such as inspiralling compact binaries. To pre- 
pare these experiments the required theoretical work consists of (i) carrying 
out a sufficiently general solution of the Einstein field equations, valid for 
a large class of matter systems, (ii) applying the latter solution to specific 
systems like binaries. 

For general isolated sources the basic problem is to relate the asymptotic 
gravitational-wave form, at the location of a detector in the source's wave zone, 
to the stress-energy tensor T""^ of the matter fields. In this article we shall 
explore this problem using the post-Newtonian expansion, when the speed of 
light c +00, as the basic approximation. This approximation has proved 
in the past to constitute an invaluable tool for understanding many aspects 
of thejjirQblcms of radiation and motion in general relativity (for reviews see 
Refs. Illam.y). 

Concerning the application to compact binaries, the point is that a model 
made of two structureless point-particles, characterized by two mass param- 
eters mi and m2 (and possibly two spins), should be sufhcient to describe 
the "inspiral" phase, driven by gravitational radiation reaction, which pre- 
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cedes the final binary coalescence. Most of the non-gravitational effects, usu- 
ally plaguing the dynamics of binary star systems, are dominated by (post- 
Newtonian) gravitational effects. However the real justification of the model 
of point particles is that the effects due to the finite size of the compact bodies 
are small |^. Inspiralling compact binaries are ideally suited for applicatiorLJif 
a post-Newtonian wave generation formalism, because it has been shown uu 
that the post-Newtonian precision required to implement successively the op- 
timal filtering technique in the LIGO/ VIRGO detectors corresponds grossly, 
in the case of neutron-star binaries, to the 3PN approximation, namely 
beyond the quadrupole moment formalism. 

In this article we review (Part A) and apply (Part B) a wave-generation 
formalism based on systematic multipole decqmpositions of Einstein's the- 
ory, ancLinitiated by Blanchet and Damour following preceding work by 
Bonnor l3 and Thorne E3. The basic multipole moments in this approach are 
some appropriately defined source moments, which are some functions of time, 
supposed to describe the actual physical source. They are iterated by means 
of a post-Minkowskian expansion of the vacuum field equations (valid in the 
source's exterior). Technically, the post-Minkowskian approximation scheme 
is greatly simplified by the assumption of multipolar expansion, because one 
can consider separately the iteration of the different multipole pieces com- 
posing the exterior field. In our formalism, the radiation field at future null 
infinity is described by the standard radiative multipole moments Ii3, which 
are obtai 
moments 



s_some explicit non-linear functionals of the more basic source 
ta. In the current situation, the closed-form expressions of the 
source multipole rnonafiuts have been established in the case where the source 
is post-Newtonian tSHlJ (existence of a small parameter which can be regarded 
as a slow motion estimate ~ v/c). For post- Newtonian sources the domain of 
validity of the post-Newtonian expansion (viz the near zone) overlaps the ex- 
terior weak-field region, so that there exists an intermediate zone in which the 
post-Newtonian and multipolar expansions can be matched together ^. Fur- 
thermore, the complete expression of the gravitational field inside the source, 
valid formally up to_any post-Newtonian order, has also been determined by 
the latter matching E2l. 

The post-Newtonian wave-generation formalism permitted to compute 
the total energy flux, £, emitted by inspiralling compact binaries. C was 



"For instance the influence of tfie Newtonian quadrupole moments, Qi and Q2, induced by 
tidal interaction between two compact objects (for which ai ~ ^^'^ '^2 — is, 

by a simple Newtonian argument, comparable (formally) with gravitational effects at the 
5PN order. |— 1|— ■ 
''This is an application of the method of matched asymptotic expansions tjo. 
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completed first to the 2PN order by Blanchet, Damour and Tver El, and, inde- 
pendently, by Will and Wiseman using their own formalism E3i23. Higher-order 
tail effects at the 2.5PN and 3.5PN orders, as well as a crucial contribution of 
tails generated by the tails themselves (the sp=jCalled "tails of tails" ) arising at 
the 3PN order, were subsequently obtained Ej. Notice that unlike the 1.5PN, 
2.5PN and 3.5PN orders which are entirely composed of tail terms, the 3PN 
order involves also, besides the tails of tails, many non-tail contributions com- 
ing from the relativistic corrections in the binary's multipole_moments. These 
have been recently computed by Blanchet, Iyer and Joguet EJ. 

The dynamics of compact binaries, i.e. their equations of motion, is also 
needed in this theoretical work. The IPN approximation in the equations of 
motion is_sta,ndard cB. The 2PN termSphave been derived by Damour aiid 



Deruelle E3cZl, Kopeikin andjGrishchuk E20, Blanchet, Faye and Ponsot Ej, 



Itoh et aZ eJ, Pati and Will E3. The 3PN order hasJaeen obtained indepen- 
dently by twQpgroups : Jaranowski and Schafer Ejo, Damour, Jaranowski 
and Schafer EjO emplcpz_lJiCLjyDM-Hamiltonian formalism of geneial relativ- 
ity; Blanchet and Faye EZIH'EScj, de Andrade, Blanchet and Faye EJ compute 
directly the equations of motion (instead of a Hamiltonian) in harmonic coor- 
dinates. To the 3PN order the end results have been shown to be physically 
equivalent in the sense that there exists a unique "contact" transformation 
of the binary's dynainical variables, that changes the harmonic-coordinates 
Lagrangian of Ref. ell into a new Lagrangian whose-Legendre transform co- 
incides exactly with the Hamiltonian given in Ref. E3. This equivalence is of 
course quite satisfying. 



PART A : POST-NEWTONIAN SOURCES 



We describe the general solution of the field equations, valid for isolated 
post-Newtonian sources, in a series of mathematical statements, presented 
here without proofs. This presentation has the advantage of emphasizing some 
key results and outlining the basic "skeleton" of the formalism. We define 
h"^ — \J—g g"^ — 'n°'^ , where g"^ is the metric, ■q"'^ is the Minkowski metric, 
and g = dei{gai3)- The field equations in harmonic coordinates, dph"^ = 0, 
read 



^ |^|yo/3 ^ ^A'^^ih, dh, d^h) , (2) 
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where □ is the flat d'Alembertian operator, T"^ is the matter tensor and A"'^ 
denotes the non-hnear gravitational source term. 

2 Linearized vacuum equations : general solution 

The first result concerns the solution of the linearized field equations in the 
vacuum region exterior to the isolated source. 

Theorem 1 The most general solution of the linearized vacuum field 

equations, nh'^^ = = df^h'^^ , outside some time-like world tube enclosing 
the source, and stationary in the past (when t < —T), reads 

hf = kf + + dPip'l - v"^d-,ipj . (3) 

The vector ipf parametrizes a linearized gauge transformation. The first term 
is given by n 




(6) 

where II and Jl are two sets of STF-tensorial multipole moments, functions 
of the retarded time u = t — r / c, and satisfying the conservation laws : I = 
const, /j — const, = const. 

3 Exterior vacuum solution : post-Minkowskian iteration 

Consider the so-called post-Minkowskian expansion, or formal series in powers 
of Newton's constant G, 



'^Multi-indices are denoted by L = iii2 ■ ■ ■ if, the derivative operator is a short-hand 
for di-^ ■ • ■ 9ij ; a tensor 1^ is STF if and only if for any pair of indices ip, iq € L, we have 
J - ip - i,- - = I---iq---ip--- ^iid ^ipiq^---ip---iq--- — 0- The Superscript (p) refers to the time 
derivation. 
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+ 00 

c^, = ^G"/.^^^ (7) 

n=l 

that is based on the general hnearized exterior solution h^^llL, Jl] of Theorem 
|l|. We. obtain a particular solution, at any post-Minkowskian order n, in the 
form u 



hf^J^VB-l[Af]+vf, (8) 

where □"J is the usual retarded integral and f^V a specific operation (defined 
by complex analytic continuation) of taking the finite part. The finite part 
is necessary because we are looking for the exterior solution, in the form of 
a multipole expansion which is singular at the origin r — 0, supposed to 
be located inside the compact-support source. The term -y^^ is a particular 
source-free retarded solution constructed in such a way that the harmonic 
coordinate condition {dph"^ = 0) is satisfied. 

Theorem 2 i The most general solution of the Einstein field equations, in the 
vacuum region outside the source, admitting post-Minkowskian and multipolar 
expansions (and stationary in the past), is given by the construction 
It depends on the source multipole moments Il{u) and Jl{u), and on four 
supplementary functions parametrizing the gauge vector in Eq. (Q). 



4 General structure of the exterior solution 



The solution is valid in all the (weak-field) exterior region of the source. We 
investigate its expansion when r and r +oo. 

Theorem 3 i The structure of the expansion of the post-Minkowskian exte- 
rior metric in the near-zone ( when r 0) is of the type ^ .• VA^ G N, 



hf{^,t) = E^^^"(1^0^i^r,i,p,„W +«(^^) : (9) 

where m £ Z,, with mo < m < N (and ma becoming more and more negative 
when n grows), p G N with p < n — 1. The functions F^^^ ^ ^ are multilinear 
functionals of the source multipole moments. 

"^We denote by (or sometimes n<i>) the STF projection of the product of unit vectors 
"n-L = "ii ■ ■ ■ 'T'i; ; S'li expansion into STF tensors = h]^(9, <p) is equivalent to the usual 
expansion in spherical harmonics Yi^{d, <f>). Similarly we pose = Xi-^ ' ' ' ^ii = ^'"i^l and 
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This theorem shows in particular that the general form of the post-Newtonian 
expansion (i.e. when c +00) is necessarily of the type ^ (In c)f/c™. 
Theorem 4 The most general multipolar and post-Minkowskian solution, 
stationary in the past, (i) admits radiative coordinates (T, X), for which the 
expansion at infinity, R +00 with U = T — R/c = const, takes the form : 

(X, T)^Y1 ^<kjU) + O (-i^) , (10) 

where the functions K^^f, „ are computable functional of the source multipole 
moments, (ii) is asymptotically simple in the sense of Penrose^, perturba- 
tively to any post-Minkowskian order. 

In radiative coordinates the retarded time J7 is a null coordinate at future 
null infinity. Theorem ^ permits to connect our work with the study by Bondi 
et al c3 of the asymptotic structure of radiative gravitational fields. Most 
importantly, it permits to calculate the relations between the source multipole 
moments II and Jl and the radiative ones, say Ul and Vl, which parametrize 
the asymptotic wave form. These relations involve nmmi-iiQn-linear multipole 
interactions associated with tails and related effects t^llJO. 

5 Post-Newtonian wave generation 

Suppose that there exists a solution h"^ , valid everywhere inside and outside 
a post-Newtonian source, whose multipole expansion coincides with the pre- 
viously constructed exterior field : M{h"^) = /i"^. Our basic assumption is 
that of the matching between the near-zone expansion, r ^ 0, of the exterior 
multipole expansion : i.e. A^(/i"^), and the far-zone expansion, r —^ +00, 
of the inner post-Newtonian expansion : i.e. A4{h°'^). The notation is M 
for the operation of taking the multipole expansion, and an overline for the 
post-Newtonian (re-)expansion. As shown by the matching, the near-zone 
expansion of the outer field is actually a post-Newtonian expansion, and the 
far-zone expansion of the inner field is identical to the multipole decomposi- 
tion. The matching equation reads 

7W(^=5]n^r™(lnr)^F^"i^^(i)=A^(7^), (11) 

where we have given the common structure of both types of asymptotic ex- 
pansions, which is given, after formal post-Minkowskian re-summation, by 
Theorem ^. 
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Theorem 5 Under the hypothesis of matching, Eq. (f71[j, the multipole 

expansion of the field exterior to a post- Newtonian source reads 



Mih'^^) = TV □reJ[>((A"'')] - -r E Hr^^ -/l^^"" - ^/^) 



(12) 



i=0 

where the "multipole moments^' are given by 



(13) 



Here, t"^ denotes the post-Newtonian expansion of the stress-energy pseudo- 
tensor defined by Eq. 

The function Si{z) entering the moments is ^ 

A different forpiaUsm, alternaliye to Theorem || has been developed by Will 
and Wiseman E3, and proved u a to be strictly equivalent to ours. 

With Theorem || in hands, one can obtain the expressions of the source 
multipole moments II and Jl parametrizing the exterior field of Sections |^ 
^. These are the ones to be inserted into the linearized metric (||)-(||) and 
iterated by means of the post-Minkowskian algorithm (0)-(H). 



^(1) 



+ c^{i +i){i+ 2){2i + 5) '^'+^^'^-^^? " - "'y'/'^) ' (15) 

Jl{u)^TV J (fiy J dzeab<ii^SiyL-i>a'Sb 

' 2(1 i^iuo^; I ^^ ^i+^yL-i>ac^\!'^\{y,u- z\y\/c) , (16) 



d^{l + 2){2l + 3) 
where we have posed 

'^This function approaches the Dirac delta-function in the limit of large / : lim ■ 
5(2), and its integral is normalized to one : J^^-^ dz5i{z) = 1. 
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6 Post-Newtonian field in the near zone 

The post-Newtonian metric can be systematically computed by means of a 
specific variant of the operator of the "symmetric" potentials, 



:V^ = E>"^-'-^- (18) 



fe=0 



Theorem 6 The expression of the post- Newtonian field inside the source, 
satisfying correct boundary conditions at infinity (no incoming radiation), fol- 
lows from the matching equation alll) as 



-9l ' \ ' ' ' ■ (19) 



4G+2S(_y^^ \Af{t-rlc)-At\t^rlO 



c'^ ^ V. ] 2r 

1=0 



The "anti-symmetric" waves are parametrized by A^^^ (u) — (u) -\-TZ'^^ (u) , 
where is given by Eq. ^T^, and where ^ 

d^yVLj^ dz^iiz)MiT"^){y,u-z\y\/c) . (20) 

The anti-symmetric waves are associated with radiation-reaction effects. For 
instance the JF^'^-contribution yields the dominant radiation reaction force at 
the 2.5PN order 00, while the 7?.^''-contribution contains the grasotational- 
wave tails in the radiation reaction, which arise at the 4PN order cJ'ca. 

PART B : COMPACT BINARIES 

We recall that the forthcoming detection and analysis of the gravitational 
waves emitted by inspiralling compact binaries should necessitate the prior 
knowledge of tha^fiquations of motion and gravitational radiation up to the 
high 3PN order M. 



^We denote 7;(z) = —25i(z). The normalization of -fi{z) : J^°° dz 'yi{z) = 1 (where the 
integral is computed by analytic continuation in / g C), is consistent with that of Si(z). 
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7 The problem of self-field regularization 

It makes sense in post-Newtonian approximations to model the compact ob- 
jects by point-like particles, but we then face the problem of the regularization 
of the infinite self-field of the particles. The Hadamard regularization has been 
extensively used in this field, but it turned out, as we now review, to yield an 
unknow coefficient at the SPM-cmder (and not before that order). 

Jaranowski and Schafer Ejo noticed that, in their computation of the 
equations of motion within the ADM-Hamiltonian formulation of general rel- 
ativity, the "standard" Hadamard regularization leads to some ambiguous re- 
sults for the computation of certain integrals at the 3PN order. They showed 
that there are two and only two types of ambiguous terms in the 3PN Hamilto- 
nian, which they parametrized by two unknown numerical coefficients Ci^static 
and Wkinotic- rnn 

Motivated by this result, Blanchet and Faye L30 introduced an "im- 
proved" Hadamard regularization, based on a theory of pseudo-functions and 
generalized distributional derivatives, which is free of ambiguities, and leads 
to unique results for the computation of all the integrals occuring at the 3PN 
order. Unfortunately, this rejgnlarization turned out to be in a sense incom- 
plete, because it was found e3'E3 that the 3PN equations of motion involve 
one and only one unknown numerical constant, called A, which cannot be de- 
termined within the method.. The comparison of this result with the work of 
Jaranowski and Schafer on the basis of the computation of the invariant 
energy of binaries moving on circular orbits, revealed EH that 



41 

t^kinotic = ^ , (21) 
^static - ^ ■ (22) 

The ambiguity ^kinetic is therefore fixed, while A is equivalent to the other 
ambiguity ecstatic- r-, 

Damour, Jaranowski and Schafer recovered the value of wj;i„otic given 
in Eq. ( pi] ) by proving that this value is the unique one for which the global 
Poincare invariance of their formalism is verified. Since the coordinate con- 
ditions associated with the ADM approach do not manifestly respect the 
Poincare symmetry, they had to prove that the Hamiltonian is compatible 
with the existence of generators for the Poincare algebra. By contrast, the 
harmonic-coordinates condition preserves the Poincare invariance, and there- 
fore the associated equations of motion should be_Laxentz-invariant, as was 
indeed found to be the case by Blanchet and Faye EZlH, thanks to their use of 
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a Lorentz-invariant regularization c2l (hence their determination of Wkinetic) ■ 

More recently, the other parameter ^static was computed by Damour, 
Jaranowski and Schafer by means of a dimensional regularization, instead 
of an Hadamard-type one, within the ADM-Hamiltonian formalism. The 
result, which in principle fixes A according to Eq. (|2^), is 

1987 

a;,.,,. = ^ A^-^- (23) 

As Damour et al El argue, clearing up the ambiguity is made possible by 
the fact that the dimensional regularization, contrary to the Hadamard one, 
respects all the basic properties of the algebraic and differential calculus of 
ordinary functions ^ _ 

Blanchet, Iyer and Joguet Ej, in their computation of the 3PN radiation 
field of two point particles, used the (standard) Hadamard regularization and 
found necessary to introduce three additional regularization constants ^, k 
and C, playing a role analogous to the equation-of-motion constant A. Such 
unknowns come from the, computation of the 3PN binary's quadrupole mo- 
ment lij . It was found Eil that the total gravitational- wave flux, in the case 
of circular orbits, depends in fact only on a single combination of these three 
constants, 



6» = ^ 2k + C . (24) 
8 The equations of binary motion 

The equations are expressed in Newtonian-like form, by means of the coordi- 
nate positions y\, yl, and velocities vl = dy\/dt, vl = dyl/dt, of the bodies in 
the harmonic coordinate system. The 3PN acceleration of one of the bodies, 
say 1, reads 

CtiL C- C- \ O / 

The first term is the famous Newtonian law, 

= ~^n\, , (26) 



^For instance, the Hadamard regularization unavoidably violates at some stage the Leibniz 
rule for the differentiation of a product. 
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(where ri2 — |yi — y2| and ni2 = (yi — y2)/f'i2)- All the other terms represent 
the successive post-Newtonian corrections. The 3PN equations of motion, that 
we show below, possess the correct perturbative limit given by the geodesies 
of the Schwarzschild metric, stay invariant when we perform a global Lorentz 
transformation, and are derivable from a Lagrangian or Hamiltonian formu- 
lation (when the radiatioiiireaction term is neglected). The IPN term is the 
Einstein-Infeld-Hoffmann acceleration : 



^IPN — 



5 5 V 4- 



' 12 



Gm-i 



• 12 



12 



4(ni2Wi) - 3(^12^2) 



-{n\2V-if' — + 4(uiU2) — 2w. 



'12 



(27) 



The parenthesis indicate the usual Euclidean scalar prod 



ni2.vi, and v\ 



At the 2PN order, the result 



isiiifs^ss^ 



^2PN — 



57 G^mf TO2 69 G'^mxm\ 



' 12 



' 12 



12 



Gr7i2 



^("12^2)^ + ^(n-i2W2)^w? 



6(ni2W2)^(wiW2) - 2(wit;2)^ + T.k'^^i^^f ^2 + '^(pxV2)v\ - 2v\ 



G^mim2 /39 



' 12 



17 



("-12W2)^ 



-(ni2Wi)^ - 39(ni2Wi)(ni2W2) 



15 2 5 

-T^l ~ 7:(^1'"2) 



G2to2 



G^m^ 



G^mim2 



2(ni2Wl)^ - 4(ni2i;i)(ni2i;2) - 6(^12^2)^ - ^{viV2) + 4w. 



2(ni2i;i) - 2(ni2t;2^ 



63, X 55. X 
-^\ni2Vi) + — {ni2V2) 



Gto2 



- 6(ni2Wi)(ni2W2) + -j^{ni2V2) + (?^i2W2)wi - 4(ni2Wi)(wiW2) 



-t- 4(ni2W2)(flW2) + 4(71i2Wl)w2 ~ ^{ni2V2)v. 



^12 



(28) 
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The 2.5PN term - an "odd" term 1/c^ - corresponds to the damping effect 
of radiation reaction : 



12 



j _ 4 G^mim2 

^2.5PN - 7 13 

"3 ^12 

4 G^mim2 

FinaUy, the very lengthy 3PN coefficient is given by0i 



Gmi 52 Gma 2 
- 6 h — h 3w 

7-12 3 ri2 



{ni2Vi2)n\2 

(29) 



^3PN — 



Gm2 / 35 g 15 n4 2 15 4 

— 2— I -j-^(»^12W2) ^('T.12W2j Vj^ + —[ni2V2) (V1V2) 

^12 V^" ° ^ 



3(ni2U2)^(wil'2)^ - -7r("l2^2)^w| + 77("-12«2)^W?W2 



y(ni2i;2)^w| + ^ 
12(ni2W2)^(wiW2)u2 - 2(wii;2)^W2 

15, n2 4 , N 4 6\ : G"^ 17111112 ( 111, .4 



2-("i2^;2)^«t + A{v,V2H - 2^;2^ j + ^1-^ - ^(nl2^'l)^ 
171 723 383 

— ("•12Wl)^(»^12W2) ^(«12^^l)^("-12l'2)^ + — ("-12Wl)(ni2t'2)^ 

455, ,4 , 229, ,2 2 205, , , , 2 

-3-("i2W2) +— p(»^12Wl) —[ni2Vi)[ni2V2)v^ 

o 4 z 

191 91 229 

— p("l2W2)^f? - -^Wl ;^(ni2t'l)^(wiW2) + 244(ni2Wl)(71i2W2)(^'lf2) 

4 o Z 



225, ,2. N 91 2/ N 177, ,2 229, ,2 2 

— ("12'U2) (W1W2) + yWi(WiW2) ^(W1W2) +— (ni2Wl) V2 

283. w . 2 , 259^^ _^2„,2 91^,2„,2 , ^„,2 81 4 

o ^2 



2 ("i2Wi)(ni2W2)w2 + — ("-12W2) ^2 - ^"l'"2 +43(wiW2)w^ 

— 3-^ - 6(ni2Ui)^(ni2'y2)^ + 12(ni2i;i)(ni2i'2)^ + 6(ni2W2)^ 
^12 V 

4(ni2Wi)(ni2i;2)(wii'2) + 12(ni2t;2)^(wii;2) + 4(i;iU2)^ 

4(ni2Wl)(ni2W2)w2 ^ 12(«'12W2)^l'2 ~ 8(flW2)w| + 4^2 

G^rrt' ( 43 \ 

^ - {n\2V\f + 2(71i2Wl)(ni2i;2) + —{nx2V2Y + 18(wiW2) - 9^2 

^12 V ^ / 

G3mim2/415, ,2 375, , 1113, ,2 

4 -^("l2Wl) ^(»^12Wl)(»^12W2) H ^('^12W2) 

r?, V 8 4 8 
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G^mlm2f 45887, ,0 24025, 
+ — —4. I - -^{ni2Vi) + -^-{ni2Vi){ni2V2) 

10469, ,2 48197 2 36227, , 36227 2 

-(^12^2) + „,„ fi T7;7r(^-'l^2) + Q ■„ •i^2 



42 ^ 840 ' 420 ' ' 840 

+ 110(ni2i;i2)'ln (^^) - 22vl2\n 
G'^mlm2 f 3187 44 /ri2 

C^mlml f 34763 44, 41 , , 
+ ^ - -A - -TT^ + 16 



''12 



'12 \ 210 3 16 / 
G^mimi/ 10478 44^ 41 2 44, /ri2 

Gto2 / 15 ^4 45 x5 3 3 2 

-:^(ni2Wi)(ni2W2) - -3-(ni2W2) - ;r(ni2%) ^i 

+ 6(ni2t;i)(ni2V2)^(wit^2) - 6(ni2V2)^(t^it^2) - 2(ni2V2)(t^it^2)^ 

- 12(ni2t;i)(ni2i)2)^vl + 12(ni2W2)^t^i + {ni2V2)vlvl 

- 4:{ni2Vi){viV2)v2 + 8(ni2t;2)(viV2)v| 

+ 4:{ni2Vi)v2 - 7(ni2t;2)i'| 



H 3-^ ( - 2(ni2t'i)^(ni2V2) + 8(ni2t'i)(ni2t;2)^ + 2(ni2ti2)" 

''12 V 

+ 2(ni2t;i)(t;it;2) + 4(7112^2) (V1U2) - 2{nY2,V])v\ - 4{nv2V2)v\ 

G^mim2 / 243, 565, ,0, 

H 3 ^(ni2^^i) +— r(ni2t;i) (ni2W2) 



'12 



4 V - 4 



269, ,2 95, 207, ,2 
^(ni2t;i)(ni2V2) - -7:('^i2'y2) + —^\n\2V\)v^ 

137 

-^{nx2V2)v\ - 3'Q{nx2V\){vxV2) 

o 

27 81 83 

+ ^(ni2W2)(^^1^^2) + -^{nx2V\)v\ + —{nv2V2)v\ 



' 12 



4(ni2Ui) + 5(ni2i;2) 



GR16: submitted to World Scientific on February 7, 2008 13 



G^miml / 307, , 479, ^ 123 , ,2 
4 ^(?^i2Wi) + — (ni2f2j + ^(ni2Wi2j7r 

G^mfma /31397, , 36227, 

-(?^12Wl) T7^{ni2V2) 



• 12 

44(ni2Ui2) In 



V 420 



420 



^12 



(30) 



The latter harmonic-coordinates equations depend at the 3PN order on two 
arbitrary constants r'l and r'2 parametrizing some logarithmic terms, however 
these constants are not physical in the sense that they can be removed by 
a coordinatejjxanaformation. By contrast, the equations of motion in ADM- 
coordinates aoo do not involve logarithms and are free of such type of 
constants. 



9 Energy of circular-orbit binaries 



Most inspiralling compact binaries will have been circularized at the time 
when they become visible by the LIGO and VIRGO detectors. In the case 
of circular orbits the equations (p5|)-(|30|) simplify drastically. Also we can 
translate the origin of coordinates to the binary's center of mass by imposing 
that the mass-dipole be = (notation of Part A). Mass parameters are the 
total mass m = mi + 7712, the reduced mass fi = mim2/m, and the useful 
symmetric mass ratio 



_ _ TOi77l2 

m (toi + TO2)^ ' 

such that < < ^, with v = j \i\ the case of equal masses, and ^ in 
the "test-mass" limit for one of the bodies. The energy E of the binary in 
the centepjpf-mass frame is deduced from the 3PN harmonic-coordinates La- 
grangian eII and expressed in invariant form (the same in different coordinate 
systems) . This is achieved with the help of the angular orbital frequency lo of 
the circular orbit, or, more conveniently, the parameter defined by 



(31) 



(32) 



After replacement we discover that the two gauge constants r'l and r'2 cancel 
out. The end result is 
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2 

675 

^~64 



ly' X - 

4 12 



27 



209323 
4032 



205 2 110, 

vr^ A 

96 9 



19 1 2 \ 2 

H 1/ zy^ X 

8 24 

155 2 35 

u u ; 

96 5184 



.(33) 



As an application of Eq. ( p3|) let us determine the location of the niinimmn 
of the energy i? as a function of the orbital frequency for circular orbits l3, 



say 



NUMERICAL 
* 

3PN 
• 2PN 



-0.03 - 



1 PN 



Figure 1. Results for i?niin in terms of u)^i^ for equal-mass binaries (y = -j). The asterisk 
marks the result calculated by numerical relativity. The point 3PN corresponds to the value 
for the ambiguity parameter given by Eq. (p3[). 

In Fig. |l| we display the results given by the successive IPN, 2PN and 3PN 
orders [the IPN and 2PN orders are defined by the obvious truncation of Eq. 
(p3|)], and compare them with the prediction recently calculated numerically 
by Gourgoulhon, Grandclement and Bonazzola cZlc3 under the assumption of 
exactly circular orbits (helicoidal symmetry). As we see the post- Newtonian 
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and numerical predictions agree reasonably well. The post-Newtonian approx- 
imation seems to converge well toward the "exact" solution, even at this very 
relativistic stage where the orbital velocity reaches vjc ~ {Gvulli j c?)^!'^ ~ 0.5. 
However, the IPN order is clearly not precise enough at this stage. 



10 Flux of circular-orbit binaries 



The calculation of the flux necessitates, following the formalism described in 
Part A, the expressions of the source multipole moments and J^. The 
crucial moment to be calculated is the mass-quadrupole, which is to be de- 
termined with the maximal poat-Newtonian precision. It has been obtained 
by Blanchet, Iyer and Joguet cJ, at the 3PN order in the case of circular or- 
bits 0. We compute the time-derivatives of the moments with the help of the 
3PN equations of motion (|2^)-(|3C|), replace them into the total energy flux 
developed at the 3PN order, and add ajithe non- linear tail and related effects 
found up to the 3.5PN order in Ref. Ej. We must also express the result in 
terms of the frequency-dependent parameter (^). Finally we get 



C = 



32d 
'5G 







1247 
336 


44711 


9271 


65 
^+18 


9072 


^ 504 



35 

1' I X 

12 



Attx^/^ 

8191 
~"672" ~ 



535 
'24 



V nx 



5/2 



/ 6643739519 ^ 16 
V 69854400 ^ ~ 
11497453 41 



2 1712 856 

— vr C Inllo; 

3 105 105 ^ 



176 



272160 48 
16285 176419 



A 



504 



1512 



9 

19897 
378 



94403 
3024 ' 



775 
324' 



] nx'/^ 



(34) 



We have = ^ + 2k + ( (and C ~ 0.577 is the Euler constant). The flux 
(^J) agrees perfectly, in the test-mass limit ^ 0, with the-jresult of linear 
black-holes perturbations calculated by Tagoshi and Sasaki e3. 

As an application of Eqs. (Q) and (^3|) one can obtain the orbital phase 
(f> = J udt of an inspiralling compact binary at the 3.5PN order simply by 



iG 3PN quadrupole is found to involve three ambiguity constants ^, k and ^ [see Eq. 
] , as well as some rg coming from the definition of the finite part process J^V employed 
'art A, and some t'q which is in fact the "logarithmic barycenter" of the constants r[ 
and r'2 present in the equations of motion (|3o|). Both ro and r'g disappear from the result 
at the end, in agreement with general expectations. 
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using the energy balance equation 



dE 
Itt 



(35) 



The orbital phase is the crucial observable to be implemented in the construc- 
tion of accurate templates for the detection and analysis of inspiralling-binary 
signals in the LIGO and VIRGO experiments. 
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